Simulation of magnetoresistance in disordered ultracold atomic Bose gases 
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Anderson localization was first investigated in the context of electrons in solids. One of the 
successes was in explaining the puzzle of negative magneto-resistance - as early as the 1940s it had 
been observed that electron diffusion rates in some materials can increase with the application of a 
magnetic field. Anderson localization has now been demonstrated in ultra-cold atomic gases. We 
present a theoretical study of the two-dimensional ultra-cold Bose gas in the presence of disorder, to 
which we apply a synthetic magnetic field. We demonstrate that, in the ballistic transport regime 
this leads to positive magneto-resistance and that, in the diffusive and strong localization regimes, 
can also lead to negative magneto-resistance. We propose experimental scenarios to observe these 
effects. 
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The study of disorder induced localization in ultracold 
atomic gases is now well established with strong local- 
ization observed in one-dimensional (Id) quasi-periodic 
lattices [1 and Anderson localization [2j [3] in both onej4] 
and three dimensional (3d) geometries [5] with disorder 
induced through laser speckle [6j [7]. In general, local- 
ization is always expected in a disordered Id system, 
whereas in 3d there exists a mobility edge [8 between 
localized and extended states and a quantum phase tran- 
sition between metallic and insulating phases can be ex- 
pected [9 . Two dimensions (2d), as is often the case, is 
the marginal dimension between these behaviours and, 
in the solid state, has lead to interesting debate [TO] re- 
garding the potential observation of a metalic phase in Si 
MOSFETs. The observation of Anderson localization in 
2d in ultracold gases is also complicated fTl J by the possi- 
bility of impurity potentials leading to classical trapping 
of the gas when the intensity of the speckle is sufficient to 
induce localization. Studies of localization in disordered 
ultracold 2d gases are therefore of timely interest. 

Anderson localization is a single-particle interference 
phenomenon and is strongly enhanced in 2d by the in- 
creased occurance of crossing trajectories (in a path in- 
tegral picture) over and above 3d. Crossing paths al- 
ways result in closed loops that constructively interfere 
back at the origin of the path with their time-reversed 
equivalent. This enhances the probability that a par- 
ticle that starts at point r remains at point r, i.e. is 
localized. If one introduces a magnetic field (of any ori- 
entation in 3d, or with some component perpendicular 
to the plane in 2d) then this time-reversal symmetry is 
broken and the enhancement of localization destroyed. 
This is the origin of negative magneto-resistance, which 
was a thirty year puzzle until explained in the context 
of Anderson localization [T2] in 1980. The observation 
of the analogue of such negative magneto-resistance in 
an ultracold atomic gas localized by disorder would be 
unambiguous evidence that the localization was an inter- 



ference phenomenon and not classical trapping or inter- 
action induced self-trapping. 

Of course, ultracold atoms are charge neutral, so we 
cannot simply impose an external magnetic field to break 
the time-reversal sysmmetry. We can however introduce 
a synthetic magnetic field by rapidly rotating the system 
[T3] or, more practically, through the use of spatially de- 
pendent light fields to couple between internal states of 
the atoms [T4HT6] . 

In this Letter we examine a 2d ultracold Bose gas in 
an optical lattice with quasi-periodic disorder induced by 
a weak second lattice of incommensurate wavelength to 
the first. We then impose a synthetic gauge field based 
upon the Raman scheme [17] to simulate an applied mag- 
netic field, breaking time-reversal symmetry. We demon- 
strate that, in the ballistic transport regime this leads to 
(normal, positive) magneto-resistance and in the diffusive 
(weak localization) and insulating (strong localization) 
regimes can induce negative magneto-resistance. 

Our system is well described by the Bose-Hubbard 
Hamiltonian which in the absence of inter-particle inter- 
actions takes the form: 

H ^ ^ Jn^m&n&rn H~ ^ ^ ^n^n^n (I) 

(n,m) n 

where a n is the Bose annihilation operator for the n th 
site, J n ,m represents hopping from the n th to the m th 
site and (n, m) indicates the nearest neighbours m of the 
n th site[T8j[19]. 

Disorder is introduced via interference with a weak op- 
tical lattice that is incommensurate with the primary. 
This, along with an external harmonic trap, is included 
in the on-site energy term: 

e n = Vdis [cos(47rx n /A 2 + <j> x ) + cos(47n/ n /A 2 + <j>y)] 

+^rap(*n+2/n)Ml, 
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where Vdi s represents the strength of the secondary lat- 
tice, Vtrap the strength of the harmonic confinement and 
Ai (A2) is the wavelength of the primary (secondary) lat- 
tice. 

To include the effect of the synthetic gauge field we 
first assume that there still exists an orthonormal set 
of basis states that are localized to each site (modified 
Wannier functions), then impose gauge invariance for any 
observable on the Hamiltonian[20-22 . This motivates 
the modification to the hopping term in 0: 



Jn,m ^ Jfi rn Jn,m CXp 



iq 
' h 
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where the integral is taken over the shortest path sepa- 
rating the end points. With this it can be shown that any 
gauge transformation with a function f(x,y), such that 
A — >• A' = A + V/, results in a n — >• a' n = exp [—iqf/h] a n 
and conserves the site number operator (n' n 
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n n ), therefore satisfying gauge invariance. 



The spatial dependence of the hopping terms can be 
calculated analytically, if one approximates the Wannier 
functions of the tight binding model [23] with Gaussians, 
or numerically with the (appropriately modifed, in the 
presence of a field) Wannier functions themselves. In the 
absence of a magnetic field it is our experience [24], and is 
well established in the literature [25j [26], that the spatial 
dependence of the hopping term is of negligible impor- 
tance compared to the on-site energy fluctuations. We 
therefore assume that this remains true in the presence 
of a magnetic field and set J U)m = J. 

We then apply the Mean-Field approximation (a n ~ 
(a n ) = z n : z n G C) to the Heisenberg equation of mo- 
tion for the Bose annihilation operator. It is convenient 
to relabel the sites according to their spatial position, 
we therefore define z^h to be the amplitude of the site 
located at (x,y) = (k,j)ai : j, k G Z, a\ = Ai/2. The 
discrete mean-field equations of motion in the Landau 
gauge (A = —Byx.) are then given by: 



6 %j,k — l 6 ^ Zj^ k -\-l -\- Cj^Zj^fc 



with: 



= A [cos(27raj + <fr y ) + cos(27ra/c • 



+^trap(j 2 +^ 2 ), 



(3) 
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where the dot indicates the derivative with respect to 
scaled time r = t/to, to = J/h. A = Vdi s /J is the sec- 
ondary lattice parameter and v t r&p = Vtr&p/J the har- 
monic trap parameter, cp is the phase accumulated in 
the Aharanov-Bohm effect when a charged particle cir- 
cumnavigates one lattice plaquette in the anti-clockwise 
direction. 

If we neglect the synthetic gauge field this system is 
separable into two Id systems, each of which, if we set 



the lattice wavelength ratio a to an irrational number, 
is equivalent to the Aubry- Andre model [27]. Therefore, 
for A < 2 the eigenstates are extended and the dynamics 
are ballistic; for A = 2 the dynamics are diffusive; and 
for A > 2 the eigenstates are exponentially localized and 
we enter the Anderson Localized regime where, except 
for virtual transmission on the order of the localization 
length, conduction is entirely suppressed [28]. In the pres- 
ence of the synthetic gauge field, however, the equations 
are not separable and we find some interesting dynamics. 

For all simulations we initialize the wave-packet in the 
ground state of the primary lattice and a fairly strong 
harmonic trap (v trap = 10 -2 ). This ground-state has 
radial symmetry in a nearly Gaussian distribution [Fig. 
[l] (a)]. The harmonic trap is then switched off at the 
same time as the synthetic gauge field and secondary 
lattice potential are switched on. The changing mag- 
netic field gives the wave-packet a gauge dependent mo- 
mentum kick[29 . In the symmetric gauge the wave- 
function remains unchanged as the contributions from 
the two dimensions cancel. We work in the Landau 
gauge, however, so we must transform the wave- function 
using zf ik = e i(pjk / 2 z? k . 

The periodicity of the primary lattice introduces a 
cosine dispersion relation in the first Brillouin zone. 
The lowest momentum states fall on the approximately 
quadratic part of the dispersion relation and so the mo- 
tion of an initially low energy wave-packet mimics that 
of a wave-packet in free space but with a reduced effec- 
tive mass [30 . The motion is termed ballistic as the rms. 
half-width scales as t with time. 

In this picture, we can understand our system in the 
presence of a synthetic gauge field. A classical particle in 
a magnetic field would complete closed circular loops in 
the 2D plane. In our simulations, an initially Gaussian 
distribution expands to a larger Gaussian with a ragged 
boundary that appears to be rotating. The size to which 
the Gaussian can expand is determined by the starting 
size and the size of the classical cyclotron orbits. So 
far, this system is analogous to that described by [31] 
which results in the celebrated Hofstadter butterfly, also 
investigated recently [32] in the context of atomic gases. 

In the extended regime (A < 2), the wave-packet re- 
leased in to a negligibly weak magnetic field gains the 
square symmetry [33] of the reciprocal lattice, with the 
real space density at long times reflecting the initial mo- 
mentum distribution, [Fig. [I] (b)] and the rms. half- 
width follows the expected t dependence [Fig. [3]- trian- 
gle] characteristic of ballistic expansion. 

When the magnetic field is stronger, time-reversal sym- 
metry is broken and hence the a x and a y symmetries [34 
(reflection in the x — and y = planes respectively) 
are also broken. The C2 = cr x a y symmetry (rotation by 
7r) remains, however, and so the solid points and unfilled 
squares in Fig. [2] are each self-similar and the Bravais 
lattice and hence the reciprocal lattice is now rotated by 
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FIG. 1: Site occupancy (logarithmic) of (a): the initial wave- 
packet at t = and (b-c): after propagating for 1000£os in 
the secondary lattice, with A = 1, and synthetic gauge field 
(b) ip = 0.001; (c) ip = 0.1; (d) ip = 1. 





• □ • □ 

•-r--ff-:;- 

□ • □ • 

v t 

• □ | • □ 



x 

FIG. 2: Diagram of the primary lattice depicting the phase 
accumulated around a lattice plaquette in the synthetic gauge 
field and the symmetries, a x and <j y , which are broken. Filled 
and unfilled dots indicate the loss of self-similarity between 
these sets of sites due to symmetry breaking. 

7r/2. The lattice parameter is also increased by a factor 
of y/2 and hence the 'volume' of the Brillouin zone is re- 
duced by half. This results in the diamond symmetry of 
Fig. [l] (c) , a signature that would be clearly observable 
in time-of- flight type experiments. 

In a very strong gauge field [Fig. [I] (d)], the initial 
momentum distribution no longer fills the entire Brillouin 
zone and so the expanded wave-packet maintains some of 
the radial symmetry, losing the diamond-like structure. 

At long times, increasing the magnetic field has the 
effect of decreasing the wave-packet's rate of expansion 
[Fig. |3j. The classical analogy is that the particles 
are completing more of their tighter circular orbits, and 




FIG. 3: Rms. half- width of the wave-packet vs. time during 
evolution in a secondary lattice with strength A = 1 and 
synthetic gauge field: (triangle) ip = 0.001; (square) ip = 0.01; 
(cross) ip = 0.1; (circle) (p = 1. 

(Inset): Long term scaling (a oc t 1 ) of the width of the wave- 
packet. 

hence traversing less linear distance, before being scat- 
tered off the pseudo-random secondary lattice. This 
is therefore analogous with normal, positive magneto- 
rest istence. 

In the inset of Fig 3 we plot 7 = ^^^y at long 
times for a selection ofgauge field strengths. The data 
points are obtained using a temporal average about each 
point. The error bars reflect the variance in the data. 
The long time behaviour of the size of the wave-packet 
approaches t 1 ! 2 with very strong magnetic field, which is 
characteristic of diffusive expansion. 

To summarize, for weak disorder, the effect of the mag- 
netic field is therefore to change the transport from ballis- 
tic to diffusive expansion. Furthermore, with increasing 
field strength, the diffusion coefficient is reduced, consis- 
tent with positive magneto-resistance. 

In the Strongly Localized regime, A > 2, we observe 
the confinement of the wave-packet consistent with the 
Aubry-Andre model for a sufficiently weak gauge field 
[Fig. [4]- square] . This Anderson Localization is the result 
of the total destructive interference of multiply scattered 
matter- waves for any sites beyond the localization length. 

Considering any path that loops back to its origin, the 
same path traversed in the opposite (time-reversed) di- 
rection is of the same length and hence will return the 
same phase. Any closed path and its time-reversed part- 
ner will therefore interfere constructively. In the presence 
of a gauge field, however, the phase is displaced pro- 
portional to the flux enclosed which has opposite sign 
for each direction (the Aharonov-Bohm effect). The 
two paths will then interfere with an essentially random 
phase. When averaged over many paths, the backscat- 
tering of the matter- waves is now dramatically reduced, 
resulting in the destruction of the Anderson localization 
and a positive expansion of the wave-packet. 
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FIG. 4: Rms. half- width of the particle wave-packet vs. time 
during evolution in a secondary lattice with strength A = 3 
and synthetic gauge field: (square) ip = 0.001; (circle) ip = 
0.01; (cross) ip = 0.1; (diamond) ip = 0.5; (triangle) ip = 1. 
(Inset): Total particles vs. time when A = 3, ip = 1 and 
particles in the outer 25 sites of the 256 x 256 system are 
strongly attenuated, (dashed line) Full system, (solid line) 
31 x 31 sites at center of the system. 



For a sufficiently strong gauge field we clearly observe 
the destruction of the Anderson localization [Fig. [4]. 
Furthermore, increasing the magnetic field increases the 
rate of expansion of the wave-packet, characteristic of 
negative magneto-resistance. The transition occurs con- 
tinuously suggesting that Anderson Localization is bro- 
ken for any magnetic field strength, although it may be 
undetectable for the system run times that we can simu- 
late. 

It is interesting to note that localization is not com- 
pletely broken for all states. Some of the eigen-states re- 
main localized. In the inset of Fig. |4]we demonstrate the 
presence of localized eigen-states. For this plot we have 
included a large negative-imaginary term in the local part 
of the Hamiltonian for sites within 25 sites of the system 
edge (in a 256 x 256 site system). When released from the 
harmonic trap, any outward bound population is quickly 
attenuated before it can create edge effects which reflect 
back to affect the center. In the plot we observe that 
the decay of the total population saturates as all the ex- 
tended states are attenuated when they reach the edge 
of the simulation grid, leaving behind the surviving lo- 
calized states. This is most pronounced for the 31 x 31 
sites at the center of the system. Such behaviour could 
be observed experimentally by taking in situ absorbtion 
images and observing the change in particle number with 
time. 

In conclusion, we have used numerical simulations to 
demonstrate the interplay between synthetic gauge fields 
and Anderson localization in the Aubrey-Andre model. 
We observe firstly positive magneto resistance in the ex- 
tended regime and then negative magneto-resistance in 
the strong localization regime. We have demonstrated 



distinctive behaviours for each regime that should be 
experimentally observable through straightforward ab- 
sorption imaging techniques. Especially, the observation 
of negative magneto-resistance can only be explained in 
the context of an interference phenomenon. This would 
therefore be an unambiguous signature that localization, 
destroyed or reduced by the imposition of a magnetic 
field, had such interference as its origin, distinguish- 
ing it from classical trapping or interaction-induced self- 
trapping. This is especially important in 2D where am- 
biguity in the origin of localization in experiments still 
resides. 
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